Resonance terahertz detection in ungated two-dimensional electron gas 
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The response of an ungated two-dimensional electron gas(2DEG) to an external electromagnetic 
excitation is analyzed. The possibility of creating a single-mode resonance detector operating in the 
terahertz frequency range is demonstrated. 

PACS numbers: 72.30+q, 73.20Mf 



I. INTRODUCTION 

The plasma wave generation mechanism [1] in an 
ultrashort-channel Field Effect Transistor (FET) has at- 
tracted a considerable attention. The high-density 2DEG 
can be used for detection [2], [3] of high-frequency radi- 
ation as well. The generation [4], [5] and both the non- 
resonant [6] and resonant [7]- [9] detection in the terahertz 
(THz) frequency rangc(scc also [10]) have been demon- 
strated experimentally. 

Recently it has been shown that the current insta- 
bility and the resulting plasma wave generation similar 
to that in FET can be obtained in ungated 2D electron 
layers [11]. Moreover, the possibility of a single-mode 
THz generator based on an ungated 2D electron layer 
was discussed in Ref. [12]. In the present paper, we calcu- 
late the response of a detector based on ungated 2D elec- 
tron layer for an arbitrary current on the assumption of 
a dissipationless carrier transport. The role of the finite 
carrier scattering strength is discussed. We demonstrate 
the feasibility of a single-mode detector operating in the 
terahertz frequency range. 



II. ANALYTICAL APPROACH 

We use the hydrodynamic model of a weakly damped, 
incompressible charged-electron fluid placed in a rigid, 
neutralizing positive background. The mean free path 
associated with electron-electron collisions is assumed to 
be less than that related to scattering by phonons and 
impurities and than the device length. The behavior of 
the 2D fluid is described [11] by the Euler and the conti- 
nuity equations: 




FIG. 1: Ungated 2D electron gas detector setup configura- 
tion. The response A(7 = U 3 — Ud is caused by the incident 
radiation U a - 
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Here, n is the 2D electron density; v, electron flux ve- 
locity; m, effective mass; — e, electronic charge; and (f>, 
potential, so that the electric field is E = —d<fi/dx. 
Let us analyze the detector mode for the equivalent cir- 
cuit shown in Fig.l. The boundary conditions at the 
source(x = 0) and drain(x = I) are 



where j is the current fixed at the drain, and W is 
the 2D layer width. We further suppose that the inci- 
dent electromagnetic radiation with a wavelength sub- 
stantially longer than the device dimensions is collected 
by an antenna coupled with 2DEG device. The boundary 
condition at the source can be realized by implementing 
the antenna output between the source contact and the 
ground. Therefore, the background 2DEG density at the 
source n is kept constant. Then, the slot antenna pro- 
vides an ac voltage U ac — U a cos(ujt), and, hence the 
subsequent small variation of the 2DES density at the 
source n ac = n a cos(iot). Finally, connecting the drain 
to the power supply via an inductance provides that the 
current at the drain is maintained constant. 

In the absence of radiation, the steady state cor- 
responds to a current-carrying electron fluid with a 
constant density n = no( because of the local quasi- 
neutrality of the 2D fluid ) and constant flux velocity 
Wo = liE = j/eWno, where \i = er/m is the carrier mo- 
bility. The steady-state potential 0o is linear downstream 
the electron flow. In the absence of dissipation </>o = 0. 
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We now seek the solutions to Eqs.(l) and (2) in the 
following form 



n = n + rii(x, t), 
v = v + ni(x, £), 
cf> = (/) + (f>i(x,t), 



(4) 



where the time-dependent variations v±(x, t), n\{x, t) are 
small, compared to the electron fluid velocity v and 
2DEG density n time-averaged over the period 2ir/w . 
For a small input signal amplitude, n a , the perturbations 
t>i, rii, 4>\ are proportional to n a , while n — no, v — vq are 
proportional to n\. To the first order in n ai we obtain 



dt dx 
dni 



n 



m dx 
dvi 

dx 



v 



Vl 
T 

dni 
dx 



0. 



0. 



(5) 



Retaining time-independent terms of the second order 
in n a in Eq.(l) and Eq.(2), we find 
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Integration of Eq.(6) with respect to x yields the dc 
detector response A0 = <f){l) — 0(0): 
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where the last term corresponds to the conventional 
ohmic voltage drop. Looking for an explicit expres- 
sion for the dc detector response given by Eq.(7), we 
now solve Eq.(5) for small perturbations vi , rii , 4>i ~ 
exp(— iuit + ikx) and then obtain 



(u) — kvo)vi = k(f>\ Vi, 

m t 

(lu — kvo)ni = knoVi. 



(8) 



Neglecting the finite-size effects [13], we explore the re- 
lation 0i = — 2 | r fc e |" 1 known for infinite ungated 2D elec- 
tron layer, where e is the background dielectric constant. 
Therefore, the dispersion equation yields: 



(co — kvo) ( lu — kvo H — I = 2a\k\ , 



where a = 
quency f2 



(9) 



Introducing the dimensionless fre- 



LO 



scattering parameter 7 



and 



flux velocity V — vq / yla, we derive the dispersion rela- 
tion for plasma waves propagating downstream k + and 
upstream fc_ the current flow as follows: 



k+ = ±- 



ii±(n + q)v - Ji±2(n 



V 2 f 2 



I V 2 

(10) 

For dissipationless carrier transport 7 = and small 
current flux, Eq.(10) reproduces the dispersion relation 
k± = ±f£ (1 =F ^) found previously in Ref.[ll]. Then, 
in the absence of current, Eq.(10) corresponds to the dis- 
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pcrsion relationship k± — ±^j(l + ^) for a free plasmon 
in an infinite 2D media. 

Seeking the solution to Eq.(8) in the form m = 
C+ exp(ik+x)+C- exp(zfc_x) and, then using the bound- 
ary conditions described by Eq.(3), we immediately ob- 
tain 
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We can now find the detector response specified by 
Eq.(7). Nevertheless, we attempt to qualitatively de- 
scribe the detector response by using the results of 
Refs.[ll],[12] for the plasma wave instability in an un- 
gated 2DEG fluid. It is well known that the response 
of a resonance-like system to an external harmonic exci- 
tation is strongly correlated with eigenvalues of the en- 
ergy spectrum. As it was demonstrated in Ref.[ll], the 
steady-state electron flow is unstable without external 
excitation, i.e., at n a = 0. The ungated 2DEG layer can 
be viewed as a resonator cavity with the complex fre- 



quency spectrum Q 



N 



n' r 



N 



ifijy, where N 



1,2. 



the mode index. The N-th mode is unstable when the 
instability increment is positive, i.e., at Q'^ > 0. For un- 
damped plasma waves tl' N 3> the complex frequency 
was calculated in Ref.[12] 
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For a small current, Eqs.(12),(13) describe [11] the dis- 
creet frequency spectrum lu' n = ^J-kcl(2N — and 
mode-independent instability increment lu" = vq/1 — 
l/2r. 
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detector response must exhibit a single-mode behavior. 
In the opposite case of a strong dissipation O/y < 1, 
the response becomes non- resonant. In order to confirm 
our predictions, we analyze in more detail the detector 
response in two limiting cases: a) V 7^ 0,7 = and b) 



A. Selective detector 

We now find the detector response under the current- 
carrying condition 7^0 and zero dissipation r = 00. 
Substituting Eq.(ll) into Eq.(7), we obtain the detector 

response = g ^i(O) 2 — vi(l) 2 ^j as follows 



FIG. 2: Dimensionless frequency ^'^(mahi panel) and incre- 
ment fi^f (inset) vs. the dimensionless velocity V for dissi- 
pationless carriers for 7 = and N — 1,2,3 modes. The 
dotted line in the main panel represents the instability cutoff 
relation £7' = 1/2 V. The bold line represents the frequency 
and the increment (inset) of the single-mode excitation in the 
presence of a finite dissipation 7 = 0.19. 

Under current-carrying conditions, the complex fre- 
quency exhibits unexpected features discussed in detail in 
Ref. [12]. In Fig. 2, we plot both the frequency and incre- 
ment for three lowest modes as a function of the flux ve- 
locity for dissipationless carriers 7 = 0. The iV-th mode 
frequency exhibits a well pronounced red-shift upon an 
increase in the current. Then, each mode is character- 
ized by the respective increment, which is positive within 
a certain range of flux velocities < V < V^ ut , where 

V'cut _ 



2-V| 



7r(2JV-'i) * s ^ ne cu t°ff velocity of the iV-th mode. 
Actually, the current-driven cutoff of the instability is 
caused by vanishing of the upstream plasma wave group 
velocity duo/dk- when the relation Q! = 1/2 V is valid(see 
Eq.(10) at 7 = 0). We state that only the first mode re- 
mains unstable at high currents V > V 2 cut = 0.25. For 
the actual case of finite scattering (7^0), the instability 
increment is independent of the mode index(see Eq.(13)), 
and, therefore, all the modes are damped identically. It 
is easy to check whether AT-th mode is still present in the 
plasma wave spectrum by plotting the damping strength 
7/2 in the instability increment plot( see Fig. 2, inset). 
Above a certain critical value of dissipation "f cr > 0.33, 
the instability disappears. Then, for 0.19 < 7 < 7 cr , only 
the proper mode remains unstable. Note that, at small 
dissipation 7 < 0.19, an increase in the current must re- 
sult in a stcp-by-step suppression of, first, higher- and, 
then, lower- index modes. 

Now we can predict some essential features of the de- 
tector response. At lot = O/7 S> 1, the response is 
expected to be resonant at the proper- and higher-order 
mode frequencies given by Eq.(12). For intermediate val- 
ues of scattering 7 < 7 cr , an increase in the current must 
lead to suppression of, first, higher- and, then, lower- 
index resonator modes. For certain values of 7, V, the 
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where we introduced dimensionless wave vectors k± = 
k±l and F(fl) is the response function. At small fre- 
quencies S7 1 the detector response is described by the 
asymptote F(fl) ~ ^ 2 /4, and thus is independent of the 
current. Then, at a zero drain current, Eq.(14) describes 

the detector response F(Q) = ta ° 2 gf/2) w hich exhibits 
sharp resonances (see Fig. 3) at the fundamental frequency 
and its odd harmonics f2jy = Wtt(2N — 1), N = 1,2... As 
expected, in the opposite case of a nonzero current, the 
response function specified by Eq.(14) demonstrates the 
noticeable red-shift of the resonator peaks(see Fig. 3). At 
a fixed carrier flux velocity V, the detector response ex- 
hibits a cutoff at a certain frequency fi' = 5^7. When 
V > V$ ut = 0.25, only the first mode of the resonator 
can be excited, and, thus, the device operates as a single- 
mode detector(see heavy curve in Fig. 3). 



B. Dissipation induced damping of the detector 
response 

We are now interested in the case of a detector op- 
crating under zero-current conditions V = and finite 
dissipation 7 / 0. The dispersion relationship is simpli- 
fied as follows k± = ±-j (fi + £7). The detector response 



given by Eq.(7) yields 



F(fi) = 



1 



2 sinh' 



(^) +( l-^)sin 2 (f) 



7 2 + fi 2 



sinh 



(15) 



Without dissipation, i.e. at 7 = Eq.(15) we reproduce 
the above-mentioned result F(Q) = tan ^ which 
yields the zero-frequency asymptote as f2 2 /4. An in- 
crease in the dissipation results in a smoothing of the 
resonant peaks (see Fig. 4) which then shrink when the 
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the typical 2DEG density uq = 10 12 sm -2 , effective mass 
m = 0.042mo, and dielectric constant e = 14 of InGaAs 
lattice matched with InP. we obtain the proper mode 
frequency u)\ = 6.5THz. Then, assuming the 2DEG 
mobility /i = 10 5 V/sm s, we find the relaxation time 
r = 3.8 x 10 _12 s and the scattering parameter 7 = 0.07. 
Therefore, the proper mode is undamped since the qual- 
ity factor lot = 25 ^> 1. The device operates as a single- 
mode detector when the carrier flux velocity exceeds the 
value v = V 2 cut Val = 9.3 x 10 5 m/s. 

We now estimate the detector responsivity defined in 
Ref.[14]. 



Q. 



FIG. 3: Detector response function specified by Eq.(14) vs. 
the dimensionless frequency at V = 0.01,0.1 and 0.3(bold 
line) for dissipationless carriers 7 = 0. Thin line shows the 
asymptote at f2 <C 1 



a 




FIG. 4: Detector response function specified by Eq.(15) vs. 

the dimensionless frequency at 7 = 0.01, 0.1, 0.7, 2 for V = 0. 

Dashed line represents the, asymptote at, $7 S> 1, r ., 
actcctor quality factor becomes on tnc order of unity, 

i.e., lot = fi/7 ~ 1. In the opposite case of a strong 
dissipation and(or) long sample k"l — ^7/2 ^ 1 the 
oscillations excited at the source by the incident radia- 
tion do not reach the drain because of the damping. The 
boundary condition at the drain is irrelevant in this case 
and the response is independent of the sample length 
since F(Q) = a^ n a ~ V- At substantially higher fre- 
quencies £1 3> 7, the response becomes independent of 
the scattering strength F(fi) = 2/f2 2 . 



III. 



DISCUSSION 



Let us estimate the values of the parameters deter- 
mined in the present paper. For I = 1/im channel and 



R 



Acf> 
~SI' 



(16) 



where S = \ 2 G/4ir is the antenna aperture; G ~ 1.5, 
dipole antenna gain factor; A, electromagnetic radiation 
wavelength; / = electromagnetic field intensity; E, 
ac field amplitude of the electromagnetic wave in a vac- 
uum. Substituting Eq.(14) into Eq.(16) and estimating 
the ac voltage as U a = EX/A, we finally obtain 



R = R F(n),R 



7r 3 e 



8Gcelno ' 



(17) 



and, morc- 



where we use the relationship U a = — 2 | r fc|"° 
over, assume k = n/2l for actual first mode detection 
case. For channel length I = 1/itm, dielectric constant 
e = 14, typical 2DEG density no = 10 12 sm -2 , we obtain 
Rq = 76[V/W]. Since the peak magnitude of the response 
function ~ 10 3 , the detector responsivity is on the order 
of ~ 10 5 [V/W] which far exceeds the typical values for 
Schottky diode detectors(on the order of ~ 10 3 [V/W]). 



IV. 



CONCLUSIONS 



In conclusion, we demonstrate that a short-channel 
ungated 2D electron gas exhibits a resonance response 
to an incident electromagnetic radiation in the terahertz 
frequency range. The responsivity of the 2DEG-based 
detector exceeds that of Schottky diodes. We demon- 
strate the possibility of a single-mode resonance selective 
detector operating at terahertz frequencies. 
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